NACA TN 2850 


CO 
! — k 
OO 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


' a i 
! □ : 
cr j 

Lnl 
03 : 

al 

Ln : 


TECHNICAL NOTE 2850 


STUDY OF PRESSURE EFFECTS ON VAPORIZATION RATE 

OF DROPS IN GAS STREAMS 

By Robert D. Ingebo 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 



Washington 

January 1953 




ITS ti 'c •. V 


•' - r ?i'V 


LFL 2®1 


TECH LIBRARY KAFB, NM 


I curl LIBmAkY KAFB, NM 



0 

RATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



SBDS 


TECHNICAL NOTE 2650 

STODY OF ERESSDRE EFFECTS ON VAPORIZATION RATE OF DROPS IN GAS STREAMS 

By Robert D. Ingebo 


SUMMARY 


The determination of the rate of vaporization of a pure liquid 
from a spherical surface exposed to a gas stream of varying static 
pressure required the use of the heat-balance equation 


dm hA 
d0 “ H v 


At 


where dm/dfl is the vaporization rate, h is the heat- transfer coeffi- 
cient, A is the surface area of the drop, H v is the latent heat of 
vaporization, and At is the difference between the gas temperature 
and the surface temperature of the drop. Sensible heat transferred by 
the liquid was negligible in comparison with latent-heat requirements . 


Experimental data were obtained for four pure liquids tested over 
an air-stream static -pressure range of 450 to 1500 millimeters of mer- 
cury. Data were also obtained for methanol evaporating in constant- 
pressure streams of helium, argon, and carbon dioxide. The following 
semiempirical equation for predicting the evaporation rate of a drop 
was determined: 


dm 

ae " Hy 


2jtd 



where kg and k v are the thermal conductivities of the gas and the 
vapor, respectively, d is the drop diameter. Re and Sc are the 
Reynolds and Schmidt numbers, respectively, g is the acceleration due 
to gravity, l is the mean free path of the gas molecules, and c is 
the root-mean-square velocity of the gas molecules. This expression 
was obtained for the case of isolated drops and a very low value of 
approach-stream turbulence . 
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The evaporation-rate and the surf ace -temperature data obtained in 
this study show the heat-transfer coefficient to he independent of the 
static pressure of the system, and the effect of pressure on the vapor- 
ization rate to he determined directly hy the effect of pressure on the 
surface temperature of the drop. An equation is also presented for cal- 
culating the At value for water at a static pressure of 450 to 
1500 millimeters of mercury hy use of the At value for 1-atmosphere 
pressure and a calculated correction term At * . 


INTRODUCTION 

Numerous investigations have heen made to determine the effect of 
pressure on the vaporization rate of droplets in st ill air (refer- 
ences 1 to 4) . However, at present, vaporization-rate data, which show 
the effect of pressure on the evaporation rate of drops in streams of 
air or other gases, are unavailable . Since combustion-chamber pressures 
of l/3 to 4 atmospheres may he encountered in the operation of Jet 
engines , the effect of air-stream pressure on the vaporization rate of 
fuel droplets is of considerable importance in the study of combustible 
fuel vapor and air mixtures for jet combustors . Such mixtures Eire gen- 
erally formed by the injection of fuel droplets at a point upstream of 
the combustion zone. Thus, the concentration of fuel vapor in the mix- 
ture entering the burning zone is determined by the vaporization rate 
of the droplets. In order to determine this rate in air streams of 
varying pressure, an investigation was made at the NACA Lewis laboratory 
in which a range of pressure conditions encountered in aircraft combus- 
tion systems was studied. 

In reference 5, an equation is presented for simultaneous heat and 
mass transfer which showed the effect of molecular mass transfer and 
turbulent momentum transfer on the heat-transfer coefficient. This 
heat-balance equation may be used to predict the effect of drop diam- 
eter, fluid velocity, and temperature on vaporization rates. However, 
since the equation was determined for an air-stream pressure of 
740 millimeters of mercury, it cannot be applied directly to vaporiza- 
tion at different air-stream pressures . The present report was there- 
fore prepared as an extension of reference 5. 

The vaporization rates of pure liquids in gas streams having static 
pressures of 450 to 1500 millimeters of mercury were determined for 
single constant-diameter drops simulated by wetted porous spheres. The 
experiments were conducted in air streams under the conditions of con- 
stant fluid velocity with varying pressure and mass-flow rates, constant 
pressure with varying velocity and mass-flow rates, and constant mass- 
flow rates with varying pressure and air-stream velocity. The mass- 
flow rate of air was varied over a Reynolds number range of 1360 to 
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3400 based on drop diameter, and surface-temperature data were obtained 
in al 1 tests. A limited number of experiments were also made in streams 
of helium, argon, and carbon dioxide at constant pressure and varying 
mass-flow rates. 

An analysis of the problem of simultaneous heat and mass transfer 
resulted in the development of a semiempirical equation which was found 
to correlate the experimental vaporization data. This expression was 
obtained for the case of single drops and small values of approach- 
stream turbulence. 


A 

b g/w- 
c 
c ' 

Dv 

d 

dm/d0 

Fr 

f 

g 

H v 

h 

K 
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SYMBOLS 

The following symbols appear in this report: 
surface area of drop, sq cm 

molecular mass diffusivity evaluated at t av , g/(cm)(sec) 

root -mean- square molecular velocity evaluated at t s , cm/sec 

root-mean-square molecular velocity of inert gas in presence of 
diffusing vapor, cm/sec 

volumetric diffusion coefficient or diffusivity of vapor eval- 
uated at t av , sq cm/sec (reference 6) 

drop diameter, cm 

vaporization rate, g/sec 

Froude number, u 2 /dg 

proportionality constant 

acceleration due to gravity, 980 cm/sec ^ 

latent heat of vaporization evaluated at t s , g-cal/g 

heat- transfer coefficient, g-cal/(sec) (sq cm)(°C) 

mass-transfer coefficient, g/(sec)(sq cm) 

thermal., conductivity evaluated at t av , g-cal/(sec)(sq cm)(°C/cm) 
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l mean free molecular path evaluated at t s , cm 

l ' mean free molecular path of inert gas in presence of diffusing 

vapor, cm 

M molecular weight 

m quantity of material vaporized , g 

n number of molecules per unit volume 

Nu Nusselt group for heat transfer , hd/k 

Nu' Husselt group for mass transfer , Kd/bg^ w 

P static pressure, atm or mm Hg 

p^ log mean partial pressure of inert or nondiffusing gas, atm 

p difference between partial pressure of vapor at droplet surface 

an d partial pressure of vapor in air stream, or driving 
potential, atm 

R universal gas constant, 8. 31X10 ^ ergs/ (°K) (mole) 

Re Reynolds number based on drop diameter with dup and p eval- 

uated at tg and t av , respectively, dup/p 

ReSc new correlative group with dup and bg^ w evaluated at t g 

and t avJ respectively, dup/b g ^ w 

Sc Schmidt group based on mass diffusivity and evaluated at t av , 

V-fb g)V 

T temperature, °K 

t temperature, °C 

t av average film temperature, i- (t g + t s ), °C 

At difference between air temperature and Q surface temperature of 

drop or driving potential, tg - t g , C 

o 

At' temperature-difference correction term, C 

At* At value at 1 atm pressure, °C 
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u velocity of gas, cm/sec 

V molal volume , cu cm/g-mole 

Z rate of momentum transfer per unit area, g/(sec^)(cm) 

Z* momentum transferred per collision per unit area, g/(sec)(cm) 

9 vaporization time, sec 

p gas viscosity, poises 

p gas or vapor density, g/cu cm 

a diameter of gas molecule, cm 

cp prefix indicating function 

Subscripts : 

a air 

f film 

g gas 

i inertia force 

j gravity force 

m molecular scale 

n macroscopic scale 

s surface 

t temperature 

v vapor 

w water 

Superscript: 

n any exponent 
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APPARATUS AND PROCEDURE 

The apparatus used in determining vaporization rates under varying 
air-stream-pressure conditions is shown in figure 1. Air was supplied 
from the central laboratory system at 50 pounds per square inch gage 
and 16-per.cent relative humidity. The air was metered by a rotameter 
and the temperature was controlled over a range of 22° to 93° C by an 
electric resistance heater. Air leaving the heater was passed through 
a calming chamber contai nin g two sections of 200-mesh screen placed 
normal, to the air stream and through a nozzle into the test section. 

The inlet ducting was designed to minimize approach- stredm turbulence. 
Control valves placed upstream and downstream of the test section made 
it possible to adjust the air-stream pressure, velocity, and mass-flow 
rate to the desired values. The sphere was placed in the center of the 
test section and air temperatures were measured by a thermocouple 
mounted in the same horizontal plane as the sphere. Air-stream static 
pressures were measured from a pressure tap drilled in the wall of the 
test section at a point in the same horizontal plane as the sphere. 
Pressures investigated covered a range of 450 to 1500 mi l limeters of 
mercury. 

The surface temperature of each liquid drop for the pressure and 
temperature ranges studied was obtained with a fine thermocouple 
junction of 40-gage iron-constantan wire flush with the surface of the 
sphere. A fine thermocouple junction was also placed in the center of 
the sphere to obtain the liquid-core temperature. By maintaining the 
liquid temperature near the surface temperature, the Bensible-heat 
transferred by the liquid was negligible in comparison with the latent- 
heat of vaporization. 

The vaporization unit (fig. 2) consisted of a syringe, a hypoder- 
mic needle, a cork sphere mounted on the end of the needle, a rubber 
stopper, and a stainless -steel tube. A screw- feed mechanism was 
required for injecting liquid into the sphere at air-stream pressures 
other than atmospheric pressure. 

Vaporization rates were measured directly by weighing the unit on 
an analytical balance before and after each run. The technique consis- 
ted in: 

(1) 'Weighing the vaporization unit on an analytical balance 

(2) Inserting the stainless-steel-tube portion of the unit into 

the test section (fig. l) 

(3) Starting the vaporization period by withdrawing the tube 

partially through the hole in the test section opposite the 
syringe, thus exposing the drop to the air stream for 2 min- 
utes (fig. 2(b)) 
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(4) Ending the vaporization run by replacing the tube over the 

sphere and sealing with the rubber stopper (fig. 2(a)) 

(5) Removing the unit from the test section, reweighing the unit, 

and obtaining the vaporization rate as the quantity of liquid 
evaporated per unit time 

Vaporization-rate data for -air streams were obtained for the 
following conditions : (a) constant fluid velocity with varying pressure 

fl.nrl air -mass-flow rates, (b) constant pressure with varying velocity 
and air -mass -flow rates, and (c) constant air -mass-flow rates with 
varying pressure and fluid velocity. 

A limited number of tests were made at constant pressure and 
varying mass -flow rates in which methanol was evaporated into streams 
of helium, argon, and carbon dioxide. In each case, the 1-inch pipe 
was disconnected from the outlet of the air heater (fig. l) and connec- 
ted by means of pipe fittings and a copper tube to a pressure regulator 
mounted on a cylinder containing the compressed gas. Thus, the gas was 
passed through the calming chamber and the test section; evaporation 
rates and surf ace- temperature data were obtained by the method described 
for air streams. The velocity of the gas stream near the surface of the 
sphere was measured with a pitot tube connected to a micromanometer . 


RESULTS AND DISCUSSION 

Effect of Pressure on Nusselt Number 

The vaporization rate of a drop may be determined "by either the 
mass-balance equation: 


dm 

d0 


b g,w P- 


V 


Nu' 

Pb 


or the heat-balance equation: 


dm 

de 


hA . , , , . , Nu 

=- At = ltd At =- 
Hv ^ 


( 1 ) 


( 2 ) 


where p^ is the difference between the partial pressure of the vapor 
at the drop surface and in the air stream, Nu' is the Nusselt number 
for mass-transfer, and p^ the log mean partial pressure of the non- 
diffusing gas. Erossling (reference 7) presented the following expres- 
sion for the mss -transfer Nusselt group for an air temperature of 
20° C: 
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. .0.5. .0.33 / Q \ 

Nu' = 2 + 0.552 (Re) (Sc) (2a) 

However, in a more recent investigation (reference 5) it was found that 
this equation was not valid over an air temperature range of 30° to 
500° C, and the following expression was presented for the heat- transfer 
Nusselt group : 


(2b) 

where the product of the Reynolds and Schmidt numbers ReSc was defined 
as a single dimensionless n umb er . The Prandtl number does not appear 
in this expression since it remained approximately constant through- 
out the tests, but the thermal- conductivity ratio k a /k v was included 
to make the expression directly applicable to an air-vapor film. 

Equation (2b) was determined at an air pressure of 740 millimeters 
of mercury and therefore does not necessarily give the correct effect of 
pressure on the heat-transfer coefficient. It does, however , predict 

Nu = qp(Re) = <p(P) 

In order to test this prediction, evaporation rates and surface - 
temperature data were obtained for constant velocity air streams of 
12.2 feet per second and a static -pressure range of 598 to 1494 milli- 
meters of mercury. Data for these runs are given in table I and show 

that, although the pressure was varied by a factor of 2^ the Nusselt 
number remained relatively constant. Thus, the pressure effect on the 
Nusselt group predicted from the Reynolds number does not appear to hold 
for this case of simultaneous heat and mass transfer . 

The effect of air -mass-flow rate on the Nusselt number (for 
methanol) was investigated at different air-stream pressures, and data 
for these runs are recorded in table II and plotted in figure 3. From 
the two plots shown in figure 3, it is evident that the effect of pres- 
sure is not adequately expressed by the relation given in equation (2b) . 




0.5 


0.6 


2 + 0.303 (ReSc) 


Momentum-Transfer Groups 

The analogy between momentum transfer and heat transfer is simply 
the concept that both heat and momentum are transferred when a molecule 
collides with a surface or another molecule. Thus, the heat transferred 
to a vaporizing drop by a molecular collision may be predicted from a 
knowledge of the momentum transferred in the collision. From the kinetic 
theory of gases, the momentum- transfer rate between molecules and a sur- 
face may be obtained in the following manner. The rate of momentum 
transfer resulting from gas molecules striking a surface may be given 
as : 



2S 


NACA IN 2850 


9 


-2 

z i, m m PC 

where Z± m is the rate of molecular momentum transfer per unit area, 
p is the gas density, and c is the root-mean-square velocity of the 
molecules^ If this rate is divided by the rate of collisions at the 
surface e/l, where Z is the mean free molecular path, 

Z 'i,m 00 1 ( 3 ) 

where Z’p m is the momentum transferred per collision per unit area 
due to the force of inertia of single molecules. 

If groups of molecules strike a surface at a velocity which is a 
function of the main-stream velocity of the gas u, the rate of momentum 
transfer per unit area may be written 

z i,n ® P u2 

The rate of collisions at a cylindrical or spherical surface is propor- 
tional to u/d since the film thickness is a function of the diam- 
eter d. For the force of inertia of groups of molecules, the momen- 
tum transferred per collision per unit area becomes 

Z‘i,n CO pud (4) 

Consider the case of molecular groups transferring momentum due to 
the force of gravity, where film thickness is again a function of d. 

The rate of momentum transfer per unit area may be written 

z j, n ro PS* 

where p is the average density of the gets . For the force of gravity 
acting on groups of molecules, the momentum transferred per collision 
per unit area is: 

z 'j,n 00 PVs^ d (5) 

since the rate of collisions at the surface is proportional to ^/gd/d. 

For the case of momentum transfer due to the force of gravity 
acting on single molecules, the rate of momentum transfer per unit area 
may be written 

Z j,m m PS* 
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■where 2 is the mean free molecular path. The rate of collisions at 
the surface is proportional to Vgl/l, and the momentum transferred 
per collision per unit area "becomes 

Z, j,m O pVi^ l 

for the force of gravity acting on single molecules . 

The four cases of momentum transfer are summarized in the following 
table: 


Force 

Scale 

Rate of momentum 
transfer per unit 
area, Z 
(g/(sec 2 )(cm)) 

Collision 
rate at 
surface 

(collisions/sec) 

Momentum transfer 
per collision per 
■unit area, Z' 
(g/(sec)(cm) ) 

Inertia 

Molecular 

-2 

pc 

c/2 

pc2 

Inertia 

Macroscopic 

pu 2 

u/d 

pud 

Gravity 

Molecular 

pgl 

Viz/z 

P Vg2 2 

Gravity 

Macroscopic 

pgd 

jgd/d 

P Vgd d 


In a gas, there is a continuous transfer of momentum from the 
gas molecules to the containing surface. This may "be measured as the 
pressure of the gas, and the momentum transferred per collision per unit 
area pci is proportional to the viscosity of the gas . If a quantity 
of gas is forced past the surface, momentum proportional to pud is 
transferred to the surface. "When a ratio of pud to the momentum 
transferred without forced flow pel is used, the following Reynolds 
number relation can "be obtained: 

Re 03 dup/pc 2 

which is the ratio of macroscopic- to molecular-scale Inertia forces. 

If heat is also transferred across the film to the surface because of 
a temperature difference tg-t g , then the Reynolds number may be eval- 
uated at the gas temperature and another momentum- transfer ratio con- 
sidered. Molecules in the f ilm will transfer momentum proportional to 
(pcZ).j. f to the surface when heat is transferred across the film. 

Thus, ^he ratio of molecular -momentum transfer without heat transfer 
to that with heat transfer is (pcl)t,g/(pcl)t,f "Where t,g indicates 
gas temperature and t,f average film temperature. This ratio should 
be included with the Reynolds number in the case of simultaneous heat 
and momentum transfer. 

Consider the case of simultaneous momentum and mass transfer. 

Here the ratio pud/p^D^ appears to be important since groups of air 
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molecules transfer momentum equal to pud to the surface, and the 
diffusing- vapor molecules transfer momentum equal to P-^D-y away from 
the surface. The momentum transfer per collision between_air molecules 
and a liquid surface of unit area may he expressed as (p cZ)^ s pro- 
vided that there is no diffusion of vapor from the liquid surface. 
However, with evaporation and diffusion occurring at the surface, the 
effective mean free path between air molecules will be greater than Z 
since Z' varies inversely with the concentration of air molecules. 
Also, the root-mean-square velocity of the air molecules may be less 
than c since momentum is lost to the diffusing vapor. Thus, the 
collision frequency c'/l ' with mass transfer must be a function of the 
collision frequency c/l obtained with no mass transfer, but it is 
considerably smaller than c/l. 

On a macroscopic scale, Froude (reference 8) has shown that if the 
collision frequency u/d is small the effect of gravity becomes very 
important and the Froude. group (pud/p Vg3 d)^ is used to describe the 
momentum-transfer process. Thus, the molecular-scale group 
p .y/gZ z/pcZ which is analogous to the macroscopi£-scale Froude group, 
may be considered important to the process when c'/l ' is small as in 
the case of momentum transfer to a_vaporizing-liquid surface. The 
groups ReSc and (p /y/g l Z)^ s /(pcZ)^. s are therefore important in the 
case of simultaneous momentum, mass, and heat transfer. This dis- 
cussion is summarized in the following table: 


Process 

Momentum- transfer groups 

Momentum transfer 

(pcZ)i^g 

Momentum and heat transfer 

( P ua) t K (pci) ti(! 
(pcZ) t , g ' (pcl) t#f 

Momentum, heat, and mass transfer 

(pua) tjK (pcD tjg (pVi? U t> 6 

(p^)t,g (p^5t,s 


Application of Momentum-Transfer Groups to Heat and Mass Transfer 

From the analogy between heat and momentum transfer, the ratio of 
macroscopic- to molecular -scale heat transfer or the Nusselt number 
hd/k would be expected to be a function of not only the groups given 
in equation 2(b), ReSc and kg/k v , but also the molecular-scale 
momentum- transfer ratio at the liquid surface (p -/g Z Z)^ B /(pcZ).^ s . 
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Thus, equation (2b) may be rewritten 

hd/k g = cpjjteSc, kg/k T , (p VgT Z)t, s /(p cZ )t,s] ( 6 ) 

The effect of ReSc and kg/ky on the Nusselt group has already 
been shown (reference 5), and equation (6) may be written: 

hd/kg = f (kg/ky) 0 " 5 (ReSc) 0 " 6 (gZ/c 2 * ) 11 (7) 


where f is a proportionality constant and n an exponent. 

The data in tables I and II indicate that the exponent n is equal 
to 0.6. In order to verify this, vaporization-rate data for four pure 
liquids were obtained in which the air -mass-flow rate was held constant 
at 9.33 grams per second and the air-stream pressure varied over a range 
of 450 to 1500 millimeters of mercury. These data, recorded in 
table III and plotted in figure 4, confirm the fact that the exponent n 
is equal to 0.6. 

The dimensionless group gZ/c 2 may be written: 


SL = = 8 . 86X10 “ 7 — 

— p 2 2 

c 6RTjmcr Tncr 


( 8 ) 


7 

since g = 980 centimeters per second per second, and R = 8.31X10 

ergs per (°K) (mole). The number of molecules in a unit volume of 

gas at any pressure P or temperature T may be written 


n = 


6.02Xl0 2S P 
22,400 760 


™ = 0.966X10 19 | 


where P is in millimeters of mercury and T is in °K. By substitu- 
ting this value for n in equation (8): 

il = 9.18X10 -26 JL 
c 2 Pa 2 

Thus, the dimensionless group gZ/c 2 is a function of pressure, mole- 
cular diameter, and molecular weight of the gas. 


In the vaporization studies, air was used as the inert heat- 
transfer medium. The best available values of collision diameter were 
used here, but they must be regarded as approximate. Since nitrogen 
and oxygen have collision diameters of approximately 3.77X10 -8 and 
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3.62X10 -8 centimeters, respectively, (reference 9)> the diameter of an air 
molecule may "be taken as approximately 3.7X10-8 centimeters. Thus, for 
air 

2L = 29 _ 2.12X10 16 

c 2 (3.7X10” 8 ) 2 

In order to test the applicability of the group g2/c^ to gas 
streams other than air, vaporization tests were performed in which 
methanol was evaporated from the surface of a porous sphere into 
constant-piessure streams of argon, helium, and carbon dioxide. By 
this means, it was possible to test the effect of molecular properties 
separate from the pressure effect given in the group gl/c . Prandtl 
numbers and values of M/cr 2 for these gases and air are shown in the 
following table: 


Gas 

cr 

(cm) 

(references 9 and 10) 

Molecular 

weight 

M/a 2 

Prandtl 

number 

Air 

3.70X10" 8 

29 

2.12X10 16 

0.74 

Helium 

2.20 

4 

0.83 

.69 

Argon 

2.86 

40 

4.88 

.68 

Carbon dioxide 

3.40 

44 

3.81 

.80 


By obtaining evaporation rates in gases other than air, it was 
also possible to vary the thermal-conductivity ratio kg/k v by a fac- 
tor of approximately 10. Data for these tests are recorded in table IV 
and plotted in figure 5. Data given in table HI are also plotted in 
figure 5, and a straight line having a slope of 0.5, which was obtained 
in reference 5, is drawn in the figure for comparison. Thus, by use of 
the molecular-scale momentum-transfer group gl/c*", it was possible to 
correlate vaporization-rate data for methanol evaporating in streams of 
air, helium, argon, and carbon dioxide with the thermal-conductivity 
ratio for a given vapor-gaB system. 

Equation (2b) predicts that the Nusselt number for a still gas may 
be expressed as: 

, / n 0 - 5 

Nu = 2(k g /k v ) 

The data used in deriving equation (2b) (reference 5) may be replotted 
as shown in figure 6. The equation for this plot is 

Nu = 2 + 0.303 (ReSc)°' 6 (k a /k v )°' 5 


(2c) 



14 


MCA TN 2850 


This equation shows the Nusselt number for still air to be independent 
of the thermal-conductivity ratio. 

Thus, the data in reference 5 cannot be used to determine whether 
equation (2b) or (2c) is the correct form since the Reynolds numbers 
were relatively high. However, in the present investigation, evapor- 
ation data for the helium-methanol system extended the range of the 
data to include high values of the thermal-conductivity ratio at low 
Reynolds numbers, and definitely show that equation (2c) is the correct 
form. The data of tables III and IV, plotted in figure 7, yield the 
following expression for Nusselt number 


^ = 2 + 2 .58X10 6 
k g 





0 ) 


Data by Frossling (reference 7) are included in this plot to exrtend 
the correlation to low ReSc values. In reference 5, the value of 
(gZ/c^)*"*'^ is 1.125X10 -7 . Thus, the proportionality constant f is 
0.303 divided by 1.125X10 -7 or 2.69X10 6 which is approximately 5 per- 
cent higher than the value obtained for f in this investigation. 

If the terms in equation (9) for the Nusselt number are substitu- 
ted in equation (2), the following vaporization- rate equation is 
obtained: 


dm 

d0 


k§At 

Hv 


2itd 



( 10 ) 


where At is the difference between the gas and the surface temper- 
ature in °C, and H v the latent-heat of vaporization in gram-calories 
per gram. In this equation, the group (ReSc gZ/c^) is independent 
of pressure, and the effect of gas-stream pressure on the evaporation 
rate may be determined directly from the effect of pressure on the 
driving potential At. 


In order to test the applicability of equation (10) to the case of 
vaporization of a fuel drop in a more humid air stream, a series of 
tests was performed in which n-octane was vaporized in an air stream 
having a relative humidity of 48 percent. The vaporization rate and 
the surface temperature of the droplet were within 2 percent of the 
values previously determined in an air stream having a relative humidity 
of 16 percent . From this study, the relative humidity of the air stream 
appeared to have very little effect on the vaporization rate of the fuel 
droplet; this agrees with the fact that the effect of water vapor on the 
value of M/a^ for the gas stream was small. 
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Effect of Pressure on Driving Potential 

Psychrometric data for water, over a pressure range of 147 to 
1523 millimeters of mercury are shown in figure 8. Values of At' were 
obtained from figure 8 at constant air temperatures above 100° C by use 
of the relation At' =0 at a total pressure of 1 atmosphere. These 
changes in At with pressure at constant temperature may be plotted 
against the pressure in atmospheres as shown in figure 9. The equa- 
tion of the straight line in the figure is 

At' = -33 log P (11) 

where At' is the temperature-difference correction in °C and P is 
static pressure in atmospheres. Thus, the At values for water at any 
pressure from 0.2 to 2.0 atmospheres may be determined, directly from the 
expression 


At = At* -33 log P (12) 

where At* is the At value at 1 atmosphere. ‘ It will be necessary 
to relate At and P for other liquids if equation (10) is to be 
used for direct calculation of evaporation rate. 

,_2 

Application of Group gZ/c to Data on Heat Transfer Inside 

Wetted -Wall Columns 

Vaporization-rate data are reported by Gilliland and Sherwood 
(reference ll) for the evaporation of liquids from the inside of wetted- 
wall cylinders into varying pressure and mass -flow-rate air streams 
The mass -transfer Wusselt number Hu' was found to be a function of 
only the Reynolds and the Schmidt numbers provided that the air-stream 
pressure is assumed to have little effect on the liquid-surface temper- 
ature or driving potential p v (see equation (l)). In their study it 
would have been very difficult to obtain surface-temperature data. 
However, in the present investigation, these data were obtained and may 
be applied to Gilliland and Sherwood's data as shown in table V. In 
this table, a portion of their data on vaporization of water (table I 
of reference 9) has been selected for examination at air-stream Reynolds 
numbers of approximately 2150. Estimated At values obtained from the 
experimental data in figure 9 are included in table V as are calculated 
values of the dimensionless groups Nu ’and gZ/c 2 , which are also plot- 
ted in figure 10. 

In the case of evaporation from inside wetted -wall columns, the 
usefulness of the group gZ/c 2 in explaining the effect of pressure 
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on simultaneous heat and mass transfer is readily apparent from fig- 
ure 10 when hoth the Reynolds and Schmidt numbers are constant. Thus, 
their data and the data in this study show the Nusselt number for simul- 
taneous heat and mass transfer to be a function of the three groups 
gZ/c^, Re, and Sc. If the Gilliland and Sherwood data were extended to 
air temperatures of 90° C or above, using water, the effect of pressure 
on the Reynolds number would be much greater than the effect of pressure 
on the driving potential At, as shown in figures, and in this cane the 
importance of using the dimensionless group gZ/c would be shown more 
clearly. 


CONCLUSIONS 

The experimental data show that the effect of gas-stream pressure 
on the evaporation rate of a pure liquid drop may be described by the 
equation 


dm k g^t 
d0 Hy 


2 itd 


1 + 1.29X10 6 



where dm/d0 is the vaporization rate - , k„ is the thermal conductivity 
of the gas, t is the difference between the air temperature and the stir- 
face temperature of the drop, Hy is the latent heat of vaporization, 
d is the drop diameter. Re is the Reynolds number. Sc is the Schmidt 
number, g is the acceleration due to gravity, Z is the mean free mole- 
cular path, c is the root-mean-square molecular velocity, and k v is 
the thermal conductivity of the vapor j ReSc is a new correlative group, 
and gZ/c^ is the molecular-scale ratio of gravitational to inertial 
forces. Thus, the Nusselt number is independent of the static pressure, 
and the effect of gas-stream pressure on the vaporization rate may be 
■determined directly from the effect of pressure on the driving force 
At. The preceding relation was evolved for the case of an isolated drop 
and a very low value of approach- stream turbulence. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 10, 1952 
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APPENDIX - CALCULATION OF DIMENSIONLESS GROUPS 
Nusselt group for n-octane (table i) . - 


Sphere diameter, d, cm 0.724 

Air temperature, t a , °C 28° 

Surface temperature, t s , °C 18 

Temperature difference, t a - t s , °C . . . 10 

Average film temperatur e , t av , °C 23 

Vaporization rate, dm/d0, g/sec 5.50X10" 4 

Latent heat of vaporization at t s , g-cal/g ’ 89 

Thermal conductivity of air at t~, k a , 

g-cal/(sec) (sq cm) (°C/em) 5.95X10“^ 


The Nusselt number may be expressed as 


Nu 


dm/de 
ird k a At 


Thus 




Nu 


(5.50X10~ 4 )(89) 

(3 . 14 ) ( 0 . 724 ) ( 5 . 95X10 -5 ) (10 ) 


Reynolds group for n-octane (table i) . 


Mass-air-flow rate, lb/hr . . . . 

Sphere diameter, d, cm 

Area of duct, sq cm ^ . 

Average film temperature, t av , C 

Viscosity of air at tf, p a , poise 


60 

0.724 

21.2 

23 

1.825X10" 4 


The mass-flcw rate per unit area up may be calculated as follows : 

” p “ (36ro°«c/wf{S ' !a) = °- 357 

Thus, the Reynolds number is 


Be = SHE = j0 .734.) (0.357 ) , _ 1413 
(1.825X10" 4 ) 
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Schmidt group for n-octane, (table III). - 


Average film temperature , t av , °C 25 

Viscosity of air at tf, p a , poise 1.835X10"^ 


The molecular mass diffusivity hg w may he calculated from the 
Gilliland expression . (reference 6) as rollows: 


_ PM,.' 


BV 

kT 


(0.0043) 


J3/2 (l\ * 1/K)°' S 

P (v a 1/3 + V /3 ) 2 


b_ „ = (0.0043) (298) 0,5 (0.00269) = 2.77X10" 4 g/(cm)(sec) 

e>^ w o 

Thus, the Schmidt number is 

Sc = JifL_ = = o.662 

b gj.w 2.77x10 

,2 

Molecular-scale momentum- transfer group gZ/c . - 
The mean free molecular path l may be expressed as 

l = l/*/2ima 2 

1 q 

where the number of molecules per unit volume n is 2.7X10 at a temper- 
ature of 0° C and a pressure of 1 atmosphere and the diameter of air 
molecule a is 3.7X10"° cm. 

Thus 


l = 6.11x10“° cm 

The root-mean-square molecular velocity c” may be expressed as 


- 3RT 

c -fw 


c 2 = (5)(8 - 51 ^) )(275) = 2.35X10 9 (cm/sec) 2 


Thus 
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Since g is 980 centimeters per second per second. 


si m (980) (6.UX1 0 : e J „ 2 . 55X10 -12 
c 2 (2.35X10; 
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TABLE I - VAPORIZATION FROM 0 . 724-CENT3METER-DIAMETER SPHERE IN 
AIR STREAM AT TEMPERATURE OF 28° C AND VELOCITY 


OF 12.2 FEET PER SECOND 


Air-flow 

rate 

(g/sec) 

P 

(mm Hg) 

At 

(°C) 

(Wde 

(g/sec) 

|i a j evaluated 
at tav 
(poises) 

(a) . 

Re 

(g-cal/g) 
(a) 

ka 

^(sec)fsq. cim)(°C/cmj} 
(a) 

Nu 

n-Octane and air 

— 

7.57 

598 

10.0 

5.50X10‘ 4 

1.825X10' 4 

1413 

88.9 

5.95X10" 5 

36.2 

9.47 

747 

8.3 

4.53 

1.830 

1762 

88.5 

5.96 

35.6 

11.83 

934 

6.7 

3.65 

1.832 

2200 

88.4 



5.98 

35.5 

Carbon tetrachloride and air 

7.57 

598 

33.6 

30.8X10" 4 

1.772X10 -4 

1452 

52.5 

5. 76X10 ~ 5 

36.7 

9.47 

747 

30.1 

27.5 

1.775 

1818 

52.3 

5.79 

36.3 

11.83 

934 

27.0 

24.6 

1.780 

2262 

52.1 

5.80 

35.9 

14.20 

1121 

25.4 

23.3 

1.785 

2710 

52.0 

5.82 

36.4 

16.58 

1307 

24.0 

22.0 

1.790 

3150 

51.9 

5.84 

35.6 

18.93 

1495 

22.0 

20.8 

1.795 

3600 

51.7 

5.86 

36.7 


CO 

H* 


Values from, reference 12 


MCA TN 2850 
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TABLE II - VAPORIZATION OF METHANOL FROM 0 . 71-CENTIMETER-DIAMETER 


SPHERE IN AIR STREAM AT TEMPERATURE OF 28° C a 



Air- mass - 
flow rate 
(g/sec) 

P 

(mm Hg) 

At 

(°C 

dm/d8 

(g/sec) 

Re 

Nu 

Nu 

/ / >0*5 

(*a/>Sr> 

ReSc 

-4 

= 1.775X10 poises} H^ = 285 

g-cal/g; kj 

3_ = 5.79X10' 

and 

k v = 3.57X10"^ g-cal/(sec) (sq. cm)(°C/cm)} 
= 1.765X10" 4 g/cm sec. 

and 


20.55 

740 

30 

1 

t 

(M 

• 

CD 

3,870 

60.7 

47.7 

3,892 

16.65 

740 

30 

7.22 

3,136 

53.1 

41.7 

3,154 

13.00 

740 

30 

6.42 

2,448 

47.2 

37.1 

2,462 

9.33 

740 

30 

5.25 

1,757 

38.6 

30.3 

1,767 

-4 

|i v = 1.760X10 poises; H v = 286 

g-cal/g; k £ 

i. = 5.74X10' 

' 5 * 
and 

k v = 3.54X10 -5 g-cal/(sec) (sq cm)(°C/cm); 
bg^ w = 1.758X10 -4 g/cm sec. 

and 


9.33 

506 

36 

7.52XHT 4 

1,772 

46.6 

36.7 

1,775 

8.58 

506 

36 

7.08 

1,630 

43.8 

34.5 

1,634 

7.82 

506 

36 

7.02 

1,485 

43.6 

34.3 

1,488 

10.10 

506 

36 

7.57 

1,920 

47.0 

36.9 

1,923 

10.72 

506 

36 

8.07 

2,100 

50.1 

39.4 

2,104 

13.00 

506 

36 

9.40 

2,470 

58.3 

45.7 

2,475 

16.65 

506 

36 

10.15 

3,180 

63.0 

49.5 

3,185 

20.55 

506 

36 

11.98 

3,900 

74.4 

58.4 

3,905 

29.20 

506 

36 

14.28 

5,550 

88.5 

69.5 

5,557 


Values for jj.^, H^, k a , and k^ from reference 12. 




TABLE TTT _ VAPT1RT7J1TTOW TTROM RPHliRRA TW ATR STRTEAM AT A MA A.A -TFT fMJ RATH OP Q 


p 

(mm Hg) 

At 

(°c) 

dm/d0 

(g/Bec) 

Hu 

2 

(cm) 

Si 

c 2 

^ReSc gZ/c 2 ^'^ 

jReSc gZ/c 2 ^^(ka/ky) 0 ’ 5 

n-Octane and 

air; t a = 

29.4 

0 C; sphere dlam. =0.78 cm; pf = 1.835X10“ 4 poises; Re = 1875; 

= 

2.53xlO y (cm/sec)^; 

Sc = 0.662; 

H v = 88.5 g-cal/g; and k a 

= 6.06X10 -5 and 

1% = 

2.20X10“^ g-cal/(sec 

)(sq cm)(°C/cm). 



795 

8.6 

5.06X10" 4 

35.5 

6.32X10 -6 

2.45X10 -12 

7.75X10" 6 

12.80X10" 6 

795 

8.6 

5.01 

36.2 

6.32 

2.45 

7.75 

12.80 

1016 

6.6 

3.34 

30.6 

4.94 

1.92 

6.70 

11.05 

1295 

5.5 

2.25 

24.7 

3.89 

1.55 

5.90 

9.73 

I5J1 

4.9 

-i r-, /-» 

l,DO 

cz, .y 

rr rx j 
«□ . 

*i nr\ 

J, . C'O 

6 . a 7 

8.70 

892 

7.2 

3.77 

31.6 

5.64 

2.18 

7.25 

11.95 

1143 

6.1 

2.71 

26.8 

4.38 

1.70 

6.24 

10.30 

1016 

6.6 

3.13 

28.7 

4.94 

1.92 

6.70 

11.05 

638 

10.5 

6.83 

39.2 

7.86 

3.06 

8.88 

14.63 

A OC 

•3ZC/*J 

11 Q 

-L-L , LJ 

O *71 

Q « / JL. 

A A c 
m . u 

in 1 R 

J-\S • J.*-/ 

*z nC 

O 4 £7*-/ 

T A 7X 

n n r\x 

_L f • 

717 

10.0 

6.16 

37.2 

7.02 

2.70 

8.25 

13.60 

Benzene 

and air; t a = 

27.8° 

C; sphere dlam. =0.78 

cm; pf = 1.787X10 -4 poises; Re = 1925; 

E® = 

2.48X10 9 (cm/sec) 2; Sc » 0.750; 

Hr = 106.5 

g-cal/g; and ka 

= 5.82X10" 5 and 

* 7=2 

35X10 g-cal/(sec)(sq cm)(°C/ 

cm) . 



798 

24.4 

12.4X10 -4, 

38.0 

6.17X10" 6 

2.44X10" 12 

8.47X10 -6 

13.33X10 -6 

798 

24 .4 

12.9 

39.5 

6.17 

2.44 

8.47 

13.33 

79fl 

24.4 

12.5 

38.2 

6.17 

2.44 

8.47 

13,33 

892 

23.2 

10.8 

34.7 

5.53 

2.19 

7.94 

12.50 

892 

23.2 

11.0 

35.4 

5.53 

2.19 

7.94 

12.50 

495 

30.0 

20.3 

50.5 

9.95 

3.95 

11.33 

17.83 

1310 

19.3 

7.5 

29.0 

3.77 

1.49 

6.32 

9.95 


Values' of pf, H v , k a > and ky from reference 12 




CO 

OJ 
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TABLE XU - Concluded. VAPORIZATION FROM SPHERES HI AIR STREAM AT A MASS-FLOW RASE 


OF 9.53 GRAMS PER SECOHD a 




p 

(mm Hg) 

At 

(°c) 

dm/d0 

(g/sec) 

Nu 

Z 

(cm) 

gl/c 2 

(ReSc gZ/c 2 ) 0 ' 6 

^ReSc gZ/c 2 ^°’ 6 x(k a /ky)°‘ 5 

Acetone 

and air; ta = 

27.8° 

C; sphere diem. =0.64 cm; \if = 1.760X10 poises; Re = 1602; 

S 2 = 

k-y - D 

2.28X10 9 (cm/sec) 2 ; Sc = 0.83; By = 135.5 g-cal/g; and k a 
2.43X10” 5 g-cal/(sec)(sq cm)(°C/cm). 

= 5.74X10" 5 and 

1021 

34.0 

ns 

30.4 

4.44X10" 6 

1. 92X10" 12 

6.98X10" 6 

10.72X10" 6 

1036 

33.8 

8.6 

29.9 

4.38 

1.89 

6.90 

10.50 

742 

37.1 

10.9 

34.5 

6.12 

2.63 

8.44 

12.98 

742 

37.1 

11.3 

35.7 

6.12 

2.63 

8.44 

12.98 

1196 

32.4 

7.6 

27.5 

3.80 

1.63 

6.34 

9.75 

1196 

32.4 

7.2 

26.1 

3.80 

1.63 

6.34 

9.75 

890 

35.4 

9.8 

32.5 

5.12 

2.19 

7.56 

11.62 

890 

35.4 

9.5 

31.5 

5.12 

2.19 

7.56 

11.62 

890 

35.4 

10.0 

32.2 

5.12 

2.19 

7.56 

11.62 

648 

38.4 

12.8 

39.1 

7.02 

3.01 

9.15 

14.06 

446 

41.6 

18.5 

52.2 

10.20 

4.39 

11.50 

17.70 

470 

41.1 

17.0 

48.5 

9.68 

4.16 

11.05 

17.00 

Acetone 

and air; = 

93.0° 

C; sphere diam. = 0.688 cm; (if = 1.945X10" 4 poises; Re = 1551; 

c 2 = 2.38X10 9 (cm/sec) 2 j Sc = 0.90; Hy = 134.2 g-cal/g; and k a 
k v = 3.00X10“^ g-cal/(sec)(sq cm)(°C/cm). 

= 6.37X10" 5 and 

798 

89.0 

29.3X10 -4 

32.1 

5.94X10 -6 

2.44X10" 12 

8. 32X10“ 6 

12.1SX10 -6 

798 

89.0 

28.5 

31.2 

5.94 

2.44 

8.32 

12.15 

1153 

84.3 

22.8 

26.6 

4.08 

1.68 

6.65 

9.70 

592 

91.6 

35.2 

37.5 

7.98 

3.29 

9.90 

14.42 

452 

93.5 

43.6 

45.5 

10.46 

4.32 

11.67 

17.00 

495 

93.0 

40.8 

42.8 

9.55 

3.95 

11.10 

16.18 

Methanol and 

air; t a = 

22. 0 C 

_4 

C; sphere diam. = 0.688 cm; \if = 1.752X10 poises; Re = 1732; 

-2 

c *=* 
k T = ; 

>.30X10 9 (cm/sec) ; £ 
5.50X10“^ g-cal/(sec' 

c = 1.00; Hy = 285 g-c 
(sq cm)(°C/cm). 

:al/g; and k a = 

5.71X10 -5 and 

1321 

22.8 

2.33X10 -4 

23.7 

3.48XL0" 6 

1.48X10 -12 

7.00X10 -6 

8.94X10" 6 

792 

28.9 

4.32 

34.6 

5.78 

2.47 

9.52 

12.17 

1011 

25.8 

3.23 

29.0 

4.53 

1.93 

8.20 

10.48 

1179 

24.5 

2.74 

25.9 

3.89 

1.66 

7.48 

9.55 

892 

27.5 

3.75 

31.5 

5.12 

2.18 

8.83 

11.28 

508 

32.0 

3.61 

44.8 


3.84 

12.40 

15.82 


Values of nf , H v , k a , and ky from reference 12 . 
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TABLE 17 - VAPOR IZATION OF METHANOL PRCM 0.71 -CENTIMETER -DIAMETER SPHERE IN GAS STREAMS OF 


VARYING MASS-FLOW RATE 8, 


p 

At 

dm/de 

u 

Ifu 

Re 

ReSc 

gl/c 2 

(ReSc gl/c 2 ) 0 ' 6 

(ReSc gj/c 2 ) 0 " 6 (kg/k v )°‘ 5 

(mm Eg^ 

(°c) 

(g/sec) 

(cm/sec) 





Methanol and 

helium; t ff 

= 26° C; 

Mg » 

1.876X10" 4 poises; Hv = 

285 g-cal/g; kg 

- 33.6X10" 5 and 

1- 

*V = 

•z a r\\f t r\ —5 — ^ n f f \ f 

vj utu./ \ tjcu ) \ oy. 

/ On /„\ . 

l:iuA 

and 

u _ T <7A\/n a**4 ~ ^ f _ \ 

L>g^ v » O. /ua jlu g/v cm A Ht:c ^ 


740 

28 

10.22X10" 4 

556 

13.88 

340 

172.5 

1.024X10 -12 

1.41S<10" 6 

4.5X10" 6 

740 

28 

10.10 

556 

13.72 

340 

172.5 

1.024 

1.418 

4.5 

740 

28 

12.23 

766 

16.62 

468 

237.0 

1.024 

1.702 

5.4 

740 

28 

12.53 

766 



17.05 

468 

237.0 

1.024 

1.702 

5.4 

Methanol and Argon; t g = 

15.5° C 

> Mg - 

2.18X10' 4 poises; Hv = 

286 g-cal/g; k^ 

■ 3.90X10“^ and 

= 

3.36X10" 5 g-cal/l 

sec)(aq 

cm)(°C/cm); 

and 

3 gjV » 1.61X10“ g/(cm)(Bec). 


740 

26.5 

3.97X10" 4 

231 

49.3 

1335 

1810 

-12 

6.05X10 

16.7X10" 6 

18.0X10" 6 

740 

26.5 

3.98 

231 

49.5 

1335 

1810 

6.05 

16.7 

18.0 

740 

26.5 

4.35 

313 

54.1 

1935 

2450 

6.05 

20.1 

21.6 

740 

26.5 

4,33 

313 

53.9 

1935 

2450 

6.05 

20.1 

21.6 

Methanol and 

^2> = 

13.0° C; 

p_ = 1.44X10 4 poises; Hy m 286 g-cal/g; kg ■ 

-5 

3.39X10 and 

k v “ 

3.36X10" 5 g-cal/( 

sec)(sq 

cm)( Q c/cm); 

and hg^ w « 1.60X10 g/(cm)(eec) 


740 

24.0 

3.73X10” 4 

426 

58.8 

3870 

3480 

4.72X10" 12 

21. 3X10" 6 

-6 

21.4X10 

740 

24.0 

2.56 

240 

41.0 

2180 

1965 

4.72 

15.2 

15.3 

740 

24.0 

2.56 

240 

40.4 

2180 

1965 

4.72 

15.2 

15.3 


Values for |i g , H v , kg, and k v from reference 12. 
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TABLE V - VAPORIZATION OF WATER FROM INSIDE WETTED-WALL CYLINDERS 


Re 

P 

(mm Eg) 

dm 

ae 

rc) 

H v 
g-cal 

k a . 

/ g-cal \ 

Nu 

ZL 

c? 



LsJi 


g 

l(sec)(sq cm)(^C/cm) , 





^secy 


• 




(a) 

(a) 

(a) 

(b) 

(c) 

(c) 

(a) 

(a) 

2120 

249 

0.105 

32 

592 

5.94X10 -5 

89.0 

7.77X10" 12 

2080 

992 

0.250 

20 

585 

6.06 

32.8 

1.95 

2160 

765 

0.317 

24 

587 

6.06 

34.7 

2.54 

2250 

112 

2.500 

40 

597 

5.88 

172.5 

17.30 

2140 

1700 

0.133 

14 

581 

6.15 

24.5 

1.14 




a Data "by Gilliland and Sherwood (reference 11 ) . 
"^Data from fig. 9. 

Reference 12, 

^Calculated . 
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figure 1, - Bcheoatio diagram of vaporization e^uinnent 
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Figure 2. - Vaporization unit. 
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Figure 3. - Vaporization of methanol from 0.71- centimeter-diameter sphere 
in (28° C) air stream at near atmospheric and reduced pressures, and 
variable air-mass -flow rate. 
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Figure 4. - Correlation of Nusselt number with molecular- 
scale momentum- transfer group. Air-mass-flow rate, 

9.33 grams per second; static pressure, 450 to 1500 milli- 
meters of mercury. 
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Figure 5. - Correlation of vaporization data with thermal-conductivity 

ratio. 
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Figure 6. - Correlation of Nusselt number 
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Figure 7 . - Correlation of Nusselt number with momentum-transfer groups 
and thermal- conductivity ratio. 
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Figure 8. - Psychrometric plot for water at several static pressures 
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Figure 9. - Correlation of t eng) erature- difference correction 

■with static pressure. 
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10. - Correlation of Nusselt number with molecular-scale momentum-transfer 
for evaporation inside vetted-vall columns. (Data from reference 9.) 




